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NONLOCAL  EFFECTS  ON  THE  CONVECTIVE  PROPERTIES  OF  THE 
ELECTROSTATIC  CURRENT  DRIVEN  ION  CYCLOTRON  INSTABILITY 

I.  INTRODUCTION 

The  current  driven  ion  cyclotron  Instability1  (CDICI)  has  been  of 

considerable  interest  to  both  space  and  laboratory  plasma  physicists  for 

more  than  two  decades.  Recently,  we  have  demonstrated  the  importance  of 

nonlocal  effects  due  to  magnetic  shear  *  (produced  by  the  field  aligned 

current)  and  the  finite  extent  of  the  current  channel4*^.  Our  past  work 

dealing  with  these  nonlocal  aspects  of  the  electrostatic  CDICI  was 

essentially  based  on  a  normal  mode  approach.  The  local  treatment  of  the 

CDICI,1*®*'  which  assumes  an  infinite  extent  of  the  width  of  the  current 

channel,  predicts  a  finite  group  velocity  of  the  CDICI  in  both  the  parallel 

and  perpendicular  directions  relative  to  the  external  magnetic  field.  For 

a  temporally  growing  mode,  the  presence  or  absence  of  a  group  velocity  in  a 

finite  sized  region  of  space  (as  opposed  to  a  point)  makes  an  instability 

convective  or  absolute;®  these  are  two  vastly  different  physical 

scenarios.  This  fact  has  been  elucidated  by  Ashour-Abdulla  and  Kennel, ^  in 

general,  and  more  specifically  for  the  CDICI  by  Ashour-Abdulla  and 

Thorne.^  In  our  earlier  work^-^  we  noted  a  significant  difference  between 

the  local  and  the  nonlocal  results  as  far  as  growth  rates  of  normal  modes 

are  concerned.  For  instabilities  like  the  CDICI  which  are  convective,  a 

normal  mode  analysis  without  regard  to  the  convective  growth  cannot  fully 

determine  the  stability  of  a  bounded  system,  and  hence  is  somewhat  limited 

in  its  application.  Thus,  the  main  objective  of  this  paper  is  to  report 

the  nonlocal  effects  due  to  magnetic  shear  and  the  finite  current  channel 

width  on  the  convective  aspects  of  the  electrostatic  CDICI. 

Manuacript  approved  April  23, 1984. 
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II.  THEORY 

In  our  treatment  we  will  follow  the  analytical  approach  as  given  In 
Bakshl  et  al . ^  We  Introduce  a  finite  current  channel  width  In  the  x 
direction.  We  take  the  electron  distribution  function  to  be  Maxwellian 
with  a  drift  velocity  parallel  to  the  magnetic  field  which  is  x  dependent: 


V.(x>  -  V° 


where  g(5)  •  exp[-  5  J  and  x  *  x  +  v  /0  .  Assuming  that  v  /8  ~  0(p  ) 

g  y  e  y  e  e 

«  Lc,  where  Lc  is  the  characteristic  length  scale  associated  with  the 
current  channel,  we  approximate  (1)  by 


Vd(x)  -  V°g[x/Lj. 


Using  (2)  and  following  the  methods  illustrated  in  Ref.  4  we  arrive  at  the 
nonlocal,  differential  equation  for  the  CDICI: 


{U7J2  ^~2 +  Qtu»  ♦(5) "  °» 

c  35  a 


where  Q  is  defined  by, 


Q  -  -  P 


2 


(5) 


1  +  Ik  XJ  +  t  +  l  rn(b)CiZ((l  -  np)c1J 


V.(5)  V.(5) 

*'  u.  -4— >• 

e  e 


A  -  ~  ^  E  r(b)  Z((l  - 


where  b  -  k2p?/2,  V'  -  31*  /3b,  t  -  T  /T  ,  p  -  8  /<i),  C.  -  u/|k  |v  , 
y  i  n  n  ie  A 

5  -  x/Lc,  Ce  -  «/|kJve,  kg  (x)  -  kz  +  ky(Lc/LsJ  J  g(5)d5,  u  -  k((x)/ky 

and  L  is  the  characteristic  length  associated  with  the  shear  in  the 
s 

magnetic  field* 

4 

The  dispersion  equation  is  given  by: 


Q'l®.  ■  0,  (7a) 

qK  V  -  (2t  +  ljuyigt-iQiu,  etjJ 1/2,  (7b) 

where  Q'  and  Q"  are  the  first  and  second  derivatives  of  Q  with  respect 

to  5  evaluated  at  5  f  which  is  defined  by  7a.  We  solve  equations  7a  and  7b 

numerically  for  the  complex  elgenf requencles  u>,  given  the  perpendicular 

wave  vector  ky,  the  parallel  wave  vector  at  the  origin  k^,  and  other 

parameters  characterizing  the  system.  For  evaluating  the  group 

velocities  v  and  v  we  take  the  derivatives  of  the  real  part  of  the 

gz  gy 

eigenfrequency  with  respect  to  k^  and  ky,  respectively,  i.e., 

tot  /Ak  and  V  -  Aw  /Ak  .  In  the  x  direction  the  wave  energy  is  confined 
r  z  gy  r  y 

to  the  vicinity  of  5^  *  Xj/Lc« 

Figure  la  shows  the  parallel  group  velocity  normalized  by  the  ion 
thermal  velocity  v^  versus  the  ratio  of  the  current  channel  size  Lc  and  the 


Here  b 


magnetic  shear  length  Lg  for  three  different  values  of  Lg. 

0.6,  t  -  0.5,  v  ■  1837  and  V°/v  -  0.28.  In  order  to  evaluate  the  value 

d  e 

2  q 

of  [k  X  J  we  have  used  the  parameters  of  Ashour-Abdulla  and  Thorne  ,  l.e.. 


n0  10  ,  B0  **  0.06  G,  and  **  2  eV. 


For  the  first 


harmonic  (k^  )  **  0.01  and  that  term  in  equation  (5)  contributes 

negligibly.  The  most  interesting  result  of  figure  la  is  that  in  the  limit 


L_  ♦  •,  V  ♦0.  One  finds  that  as 
c  gz 


is  increased  beyond 


Lg,  approaches  zero  as  (Lg/Lc)  .  Further,  when  Lc  is  made  smaller  than 

10  L  the  value  of  V  attains  its  local  value  (l.e.,  the  value  predicted 
s  gz 

by  local  theories **”)  and  stabilizes  at  that  value.  However,  when  Lc  is 

further  reduced  such  that  L.  <  p . ,  V  once  again  decreases. 

c  i  gz 

Figure  lb  is  a  plot  of  V  against  Lc/Lg  for  the  parameters  of  Figure 

gy 

la.  Here  too  we  find  a  reduction  in  V  in  the  shear  dominated  regime 

gy 

(i.e. ,  Lc  >  Lg)  much  in  the  same  fashion  as  the  temporal  growth  rate.^ 

However,  the  difference  between  the  local  and  the  nonlocal  values  of  V  is 

gy 

q 

not  as  striking  as  that  of  V .  As  noted  by  Ashour-Abdulla  and  Thorne, 

V  for  the  CDICI  near  the  dominant  normal  mode  (which  occurs  for 

gy 

b  ^  0(1))  is  extremely  small;  and  combined  with  our  results  that  V  >0 

gz 

in  the  nonlocal  regime  the  CDICI  becomes,  effectively,  an  absolute 

instability.  For  l»c  <  p^,  becomes  negative.  Figures  2a  and  2b  are 

plots  similar  to  la  and  lb  for  the  second  harmonic.  Here  we 

use  t  ■  1,  V?/v  »  0.55  and  b  ■  2.4.  The  value  of  b  for  all  the 

q  e 

calculations  corresponds  to  the  maximum  temporal  growth  rate.  He  see 

trends  similar  to  figure  la  and  lb  for  the  higher  harmonic  also. 


In  Figure  3  we  plot  the  convective  growth  rate  k.  "  y/V  normalized 

lz  JZ 

by  the  mean  ion  Larmor  radius  versus  He  see  two  distinct  regions 


joined  smoothly.  For  very  small  L  /L.  (i.e.,  L  /L0  <  10“3),  k,  p,  <  0 

C  8  C  S  ^  1Z  1 

which  Indicates  spatially  damped  waves.  When  L  /L  is  Increased  the  curves 
attain  a  plateau  which  corresponds  to  the  local  theory.  Here  we  find  a 
finite  k^p^  a°d,  hence,  a  finite  region  in  which  the  wave  amplitude  can  be 
amplified  by  one  e-fold.  For  Lc  >  Lg,  k^p^  grows  without  bound,  implying 
that  the  wave  amplitude  can  e-fold  in  an  insignificantly  small  region 
(Recall  that  the  system  is  assumed  to  be  infinite  and  uniform  in  the  y 
direction,  hence  convection  in  this  direction  is  ignorable).  Thus  the 
character  of  CDIC1  is  changed  from  convective  to  absolute. 

We  have  so  far  considered  only  one  species  of  ions  characterized  by  a 
temperature  T^.  In  order  to  make  the  system  more  realistic,  especially  for 
space  applications,  we  added  a  second  ion  species  characterized  by  a  hotter 
temperature  but  still  with  a  Maxwellian  distribution.  The  general  features 
of  the  group  velocity  behavior  is  still  the  same  as  in  a  single  ion 
species.  A  more  detailed  parametric  study  of  the  nonlocal  convective 
aspects  of  CDICI  with  a  loss  cone  distribution  for  the  hotter  ion  species, 
along  with  the  magnetospherlc  application  will  be  presented  elsewhere. 

III.  DISCUSSION 

Thus,  in  contrast  to  local  theory,  nonlocal  theory  gives  zero  group 
velocity  for  the  CDICI  along  the  average  external  magnetic  field  direction 
in  the  limit  >  Lg.  Clearly  the  ratio  Lc/Ls  is  the  important  parameter 
and  its  magnitude  must  be  well  established  before  making  definitive 
conclusions  regarding  the  electrostatic  CDICI.  Depending  on  the  value  of 
L  /L  ,  the  CDICI  can  be  classified  into  three  main  regimes: 


(i)  L  >  0.1  L  This  regime  is  purely  nonlocal  with  a  temporal 

c  ^ _ s 

growth  rate  much  reduced  from  its  local  value.  Also  the 
Instability  is  absolute  or  only  very  weakly  convective  in  the  z 
direction  which  is  the  direction  of  the  average  external 
magnetic  field  and  the  current  flow.  There  is  no  energy  flow 
in  the  x  direction,  and  energy  flow  in  the  y  direction  is 

g 

lgnorable  if  3^  =  0. 

(11)  L  <  0.01  L  In  this  regime  nonlocal  theory  reproduces  the 

c _ s 

results  of  local  theory,  and  the  instability  becomes  convective 
in  the  z  direction  with  a  convective,  as  well  as  a  temporal 
growth  rate,  equal  to  the  ones  given  by  the  local  theory. 
There  is  a  smooth  transition  between  regimes  (1)  and  (ii)  as  is 
evident  in  all  the  figures. 

(iii)  **  The  finite  channel  width  becomes  Important. 

There  is  once  again  a  reduction  in  the  temporal  growth  rate 

(filamental  quenching^)  but  the  instability  remains  convective, 

i.e.,  V  remains  non-zero.  However,  when  L  is  further 
gz  c 

reduced  so  that  L  <  p,  both  V  and  V  become  reduced  and 

c  i  gz  gy 

eventually  become  negative. 

Physically  the  reason  for  vanishing  parallel  group  velocity  can  be 
described  as  follows.  Figure  4  shows  the  angle  space  where  the  imaginary 
part  of  the  local  dispersion  relation  Dj  <  0,  and  the  growth  rate  is 
positive  for  the  CDICI.  The  extent  of  this  unstable  angular  space  Au  for 
the  CDICI  in  our  parameter  range  is  about  0.1.  In  the  nonlocal  theory  a 
wave  packet  is  formed  whose  size  is  governed  by  the  smaller  of  the  two 
scale  lengths  Lc  and  0.1  Lfl.  When  Lc  <  0.1  Lg  the  drift  velocity  Vd(x) 
varies  sufficiently  fast  as  a  function  of  x,  or  as  a  function  of  the  angle 


u  within  the  space  Au,  so  as  to  form  an  effective  Qj  curve  which  has  the 
form  of  a  steep  well.  The  position  of  the  bottom  of  this  well  is  at  u 

o 

■  k  /k  ;  thus  different  choices  of  k  center  the  well  in  different  parts  of 
z  y  z 

the  region  Au  and  the  corresponding  wave  packets  (of  size  /p^L^  )  sample 

different  regions  of  Au  which  lead  to  different  growth  rates  Y  as  well  as 

different  real  frequencies  w  .  On  the  other  hand,  when  L  >  L  ,  the 

r  c  s 

variation  of  V^x)  in  the  domain  Au  is  rather  weak  and  the  effective  Qj 

curve  is  governed  by  the  variation  of  u  rather  than  that  of  in  Eq. 

(3).  Now,  if  we  choose  different  k  ,  an  appropriate  translation  in  x 

z 

moving  away  from  the  center  of  the  slab)  is  sufficient  to  provide 

almost  the  same  invariant  curve  as  a  function  of  u*  The  resulting  wave 

packet  (of  size  forms  at  the  same  position  in  u  space  (even  though 

we  choose  different  k  )  and  has  the  same  growth  rate  r  and  real 

z 

frequency  u>  ,  independent  of  k  .  This  makes  the  group  velocity  3cj  /3k 
r  z  r  z 

vanish  in  the  limit  L  ♦  ». 

c 

IV.  CONCLUSION 

We  have  shown  in  this  paper  that  the  nonlocal  effects  due  to  the 

magnetic  shear  (produced  by  a  field  aligned  current)  and  a  finite  channel 

current  width  can  drastically  alter  the  character  of  the  electrostatic 

CDICI  by  making  the  instability  effectively  absolute  in  the  z  direction. 

The  important  parameter  turns  out  to  be  the  ratio  of  the  two  scale  lengths 

involved  in  the  problem  L  and  L  .  Depending  on  the  value  of  L  /L  one  can 

c  s  c  s 

classify  the  CDICI  in  three  regimes  as  described  earlier,  and  a  careful 
assessment  of  the  value  of  Lc/Lg  for  a  given  physical  situation  becomes 


essential  to  draw  any  definitive  conclusion  regarding  the  growth  and 
convective  characteristics  of  the  electrostatic  CDICI.  We  would  also  like 
to  note  that  the  nonlocal  effects  due  to  magnetic  shear  can  be  expected  to 
produce  corresponding  phenomena,  including  the  vanishing  of  the  parallel 
group  velocity,  in  various  other  current  driven  instabilities  as  well. 
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ie  convective  growth  rate  normalized  by  the  Ion 


A  plot  of  Imaginary  part  of  the  local  dispersion  relation  Dg(<i>,  k) 
against  k  /k.  (reproduced  from  Ref*  2)  for  various  V^/v  .  Here  b  ■ 

Z  1  d  6 


0.5,  |i  ■  3674  and  four  different  values  of  the  temperature  ratio  t: 
(a)  t  -  0.1,  (b)  t  -  0.2,  (c)  t  -  0.33,  and  (d)  t  -  1.0. 
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